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ABSTRACT 

We present a semi-analytical model for the formation and evolution of a high redshift quasar 
(QSO). We reconstruct a set of hierarchical merger histories of a 10'^ dark matter halo 
and model the evolution of the corresponding galaxy and of its central super massive black 
hole. The code GAMETE/QSOdust consistently follows: (i) the black hole assembly via both 
coalescence with other black holes and gas accretion; ( ii) the build up and star formation his- 
tory of the quasar host galaxy, driven by binary mergers and mass accretion; ( Hi) the evolution 
of gas, stars, metals in the interstellar medium (ISM), accounting for mass exchanges with the 
external medium (infall and outflow processes); (iv) dust formation in Supernova (SN) ejecta 
and in the stellar atmosphere of Asymptotic Giant Branch (AGB) stars, dust destruction by 
interstellar shocks and grain growth in molecular clouds; (v) the AGN feedback which powers 
a galactic-scale wind, self-regulating the black hole growth and eventually halting star forma- 
tion. 

We use this model to study the case of SDSS J1148-I-5251 observed at redshift 6.4. We ex- 
plore different star formation histories for the QSO host galaxy and find that Population III 
stars give a negligible contribution to the final metal and dust masses due to rapid enrichment 
of the ISM to metallicities > = [10"^ - 1Q^*]Zq in progenitor galaxies at redshifts > 10. 
If Population II/I stars form with a standard initial mass function (IMF) and with a charac- 
teristic stellar mass of md, = 0.35 M©, a final stellar mass of [1 - 5] x lO^'M© is required to 
reproduce the observed dust mass and gas metallicity of Jl 148. This is a factor 3 to 10 higher 
than the stellar mass inferred from observations and would shift the QSO closer or onto the 
stellar bulge - black hole relation observed in the local Universe; alternatively, the observed 
chemical properties can be reconciled with the inferred stellar mass assuming that Population 
II/I stars form according to a top-heavy IMF with nich = 5 Mq. We find that SNe dominate the 
early dust enrichment and that, depending on the shape of the star formation history and on 
the stellai- IMF, AGB stars contribute at redshift z < 8 - 10. Yet, a dust mass of [2 - 6] x 10*^ 
Mq estimated for Jl 148 cannot be reproduced considering only stellar sources and the final 
dust mass is dominated by grain growth in molecular clouds. This conclusion is independent 
of the stellar initial mass function and star formation history. 

Key words: Galaxies: evolution, high-redshift, ISM; quasars: general; stars: AGB and post- 
AGB, supernovae: general, ISM: dust, extinction 



1 INTRODUCTION 

The presence of large dust masses in high redshift quasars (QSOs) 
has been revealed by mm and sub-mm observations of samples 
of 5 < z < 6.4 QSOs in the Sloan Digital Sky Survey (SDSS) 
(Bertoldi et al. 2003; Priddey et al. 2003; Robson et al. 2004; Bee- 
len et al. 2006; Wang et al. 2008). The detection of dust thermal 
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emission and the inferred far-IR (FIR) luminosities of these objects 
suggest masses of warm (T < lOOK) dust of a few IO^Mq. Observa- 
tions of QSOs and galaxies at high redshift are strongly alfected by 
the presence of dust and our ability to interpret the observed prop- 
erties of high redshift sources depends on a detailed understanding 
of the star formation history and the history of metal and dust pol- 
lution in their interstellar medium (ISM). Dust properties, such as 
grains chemical composition and size distribution, are of primary 
importance in modeling the extinction curves in both the local Uni- 
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verse (Cardelli, Clayton & Mathis 1989; Pel 1992; Weingartner & 
Draine 2001) and at high redshifts (Maiolino et al. 2004; Stratta et 
al. 2007; Gallerani et al. 2010). Yet, the origin of dust at such early 
epochs is still unclear. 

Among the stellar sources, the cool and dense atmosphere of 
Asymptotic Giant Branch (AGB) stars and the expanding ejecta of 
core collapse supernova (SN) offer the most viable sites of dust 
grains condensation. A SN origin has often been advocated as the 
only possible explanation for the large amount of dust observed 
in high redshift QSOs since it is usually assumed that low and 
intermediate-mass stars {m, < 8M0) must evolve on a timescale 
comparable to the age of Universe at redshift z > 6 (about 1 Gyr, 
Morgan & Edmunds 2003; Marchenko 2006; Dwek, Galliano & 
Jones 2007) before they reach the dust producing stages. This sce- 
nario has been tested through observations of the reddened quasar 
SDSSJ1048-f46 at z = 6.2 (Maiolino et al. 2004) and of the Spec- 
tral Energy Distribution (SED) of the z = 6.29 GRB 050904 
(Stratta et al. 2007) and the z ~ 5 GRB 071025 (Periey et al. 2010). 
In these sources, the inferred dust extinction curve is different with 
respect to any of the extinction curves observed at low z, and it 
shows a very good agreement with the extinction curve predicted 
for dust formed in SN ejecta. This is an indication that the prop- 
erties of dust may evolve beyond z > 5 - 6. In a recent study, 
Gallerani et al. (2010) find that the extinction curves of a sample of 
QSOs with 4 < z < 6.4 deviate from the Small Magellanic Cloud 
(SMC) extinction curve (the one which better describe the z < 2 
quasar reddening), suggesting that production mechanisms and/or 
dust processing into the ISM must be different at high redshift. 

On the theoretical side, the evolution of dust in the early Uni- 
verse is still poorly understood. Some of the models developed so 
far are still consistent with the scenario where most of the dust in 
QSO at z ~ 6 is produced by SNe. However these estimates either 
neglect dust destruction by interstellar shocks (Maiolino et al. 2006; 
Li et al. 2008) or make extreme assumptions on SN dust conden- 
sation factors (Dwek, Galliano & Jones 2007). In addition, it has 
been shown that only [2 - 20]% of the newly condensed dust sur- 
vives the passage of the SN reverse shock and grains are reduced 
both in number and in size, depending on the surrounding ISM den- 
sity as well as on dust composition (see Bianchi & Schneider 2007; 
Nozawa et al. 2007; Bianchi et al. 2009). 

In a previous paper (Valiante et al. 2009; hereafter V09) we 
showed that AGB stars can provide a non negligible contribution 
to dust formation at high z, being able to produce dust already 
[30 - 40] Myr after the onset of the star formation, which is about 
the mean lifetime of a < 8 Mq star (see e.g. Padovani & Matteucci 
1993; Raiteri et al. 1996). In particular, for a standard stellar initial 
mass function (IMF), the characteristic timescale at which AGB 
stars dominate dust production ranges between 150 and 500 Myr, 
depending on the assumed star formation history and on the stel- 
lar initial metallicity. Hence, we concluded that these stellar dust 
sources must be taken into account when modeling dust evolu- 
tion at high redshift. In V09 we applied a simple chemical evo- 
lution model with dust to the host galaxy of the z = 6.4 SDSS 
Jl 148+5251 QSO (hereafter Jl 148) using as an input to the model 
the star formation history obtained in a numerical simulation by Li 
et al. (2007); this simulation predicts that a large stellar bulge is al- 
ready formed at z = 6.4, in agreement with the local SMBH-stellar 
mass relation, through a series of strong bursts triggered by galaxy 
mergers in the hierarchical evolution of the system. We found that 
the total mass of dust predicted at z = 6.4 is within the range of 
values inferred by observations of SDSS JIMS, with a substantial 
contribution (80%) of AGB-dust. We concluded that stellar sources 



can account for the huge dust mass observed at these early cosmic 
epochs, even when dust destruction by SN shocks is taken into ac- 
count, at the price of having a total stellar mass of ~ lO'^Mg, which 
can hardly be reconciled with the dynamical mass of ~ 5 X IO'^Mq 
inferred from CO observations within 2.5 kpc from the QSO center 
(Walter etal. 2004). 

The above critical point has been recently emphasized by 
Michalow.ski et al. (2010) who showed that using stellar masses 
approximated as the difference between the dynamical and the H2 
gas masses for 5 < z < 6.5 QSOs, stellar sources of dust are not 
efficient enough to account for the dust mass inferred in these sys- 
tems. Therefore, they suggest additional non-stellar dust formation 
mechanism; in particular, significant dust grain growth in the inter- 
stellar medium of these QSOs. Indeed, dust grains formed in stellar 
outflows can serve as growth centers for accretion of icy mantels 
in dense Molecular Clouds (MC, see e.g. Draine 1990). It is al- 
ready well known that MC-grown dust is the dominant component 
in the Milky Way (see e.g. Zhukovska et al. 2008). Given the large 
masses of molecular gas detected in CO for most of these systems, 
Draine (2009) suggested that the bulk of their dust mass is the re- 
sult of grain growth and that SN are required to produce the metals 
that compose the grains and the "seed" grain surface area on which 
more material can grow in MCs. 

With the present work we investigate the nature of dust in 
high-z QSOs host galaxies focusing in particular on J 1148 at z = 
6.4. To this aim, we have improved the semi-analytic code GAlaxy 
MERGER Tree and Evolution (GAMETE), developed by Salvadori, 
Schneider & Ferrara (2007, hereafter SSF07) to model the hier- 
archical assembly and merger history of a QSO host galaxy at 
z = 6.4, together with its central super massive black hole includ- 
ing feedback from the active galactic nucleus (AGN) in the form 
of a galactic-scale wind. The chemical evolution network has also 
been modified to describe the evolution of dust taking into account 
the specific stellar lifetimes, and the evolution of dust in the ISM 
through destruction by interstellar shocks and dust-growth in MCs. 
This code, which we will refer to as GAMETE/QSOdust, will en- 
able us to overcome the main limitations of V09 analysis, namely: 
( i) we will be able to model different possible star formation histo- 
ries for Jl 148 discussing their impact on the chemical properties of 
the QSO host and, in particular, on the final dust mass, and ( ii) we 
will be able to explore the relative importance of stellar sources of 
dust (SN and AGB) and of grain growth in MCs. At the same time, 
the model can be constrained using the properties directly observed 
or inferred from observations of Jl 148, such as the mass of the cen- 
tral black hole, the total mass of gas and stars, the mass in metals 
and dust. 

The paper is organized as follows: in section 2 we will intro- 
duce the main physical properties of the QSO JI148 and its host 
galaxy discussing the observational uncertainties. In section 3 we 
will briefly summarize the basic features of the original GAMETE 
code and present the parameters adopted to simulate the hierarchi- 
cal assembly of a QSO host halo at z = 6.4. In section 4 we will 
discuss the new features of the code which allow to follow the for- 
mation and evolution of the centtal black hole and its effects on the 
parent host galaxy through AGN-driven feedback. In section 5 we 
describe the chemical evolution network, including the evolution 
of dust in the ISM. In section 6 we will present the results and in 
section 7 we will discuss their implications. Finally, in section 8 we 
will draw the main conclusions of this work. 

In what follows we adopt a Lambda Cold Dark Matter 
(ACDM) cosmology with 0.„ = 0.24, = 0.76, a/, = 0.04, and 
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Hq = 73 km/s/Mpc. The age of the Universe at a redshift z = 6.4 is 
900 Myr. 



2 J1148: A QSO AT REDSHIFT 6.4 

In this section we briefly summarize the physical properties of this 
extremely luminous QSO and of its host galaxy which have been 
derived observationally. The full list of properties, and observa- 
tional uncertainties, is provided at the end of this section. 

2.1 Dark matter host halo mass and AGN properties 

The mass of the dark matter halo where Jl 148 is hosted, Mj,, is usu- 
ally estimated by requiring that the number density in halos with 
masses > Mh at z ~ 6 matches the observed space density of QSOs 
at the same redshift (Fan et al. 2004). Depending on the cosmo- 
logical model, this generally yields to Mh = [10*- - 10'^^]Mq. The 
space density of observed QSOs with average redshift < z >= 6 
and luminosity M1640 < -26.7 is (6.4 ± 2.4) x lO"'" Mpc"^ (Fan et 
al. 2004). For our adopted cosmological model, this would require 
a minimum halo mass of ~ 3 x lO'^M©. To compare with previ- 
ous analyses (Volonteri & Rees 2006; Valiante et al. 2009) in this 
work we generate Monte Carlo realizations of the merger histories 
adopting a value of Mh = IO'^Mq. The implications of a different 
halo mass and redshift will be addressed in a forthcoming paper, 
where a larger sample of QSOs will be investigated. 

J1148 is a very luminous QSO with absolute AB magnitude 
of the continuum in the rest-frame at 1450 A of M1450 = -27.82 
(Fan et al. 2003). Assuming a power law fyocv" with a = -0.5 
the corresponding B-band magnitud^ and luminosities are Mb = 
-28.3 and Lg = 5.2 x IO'^Lq where the effective wavelength of the 
B-band is 4400A. 

The energy requirement to power such luminosity suggests 
the presence of a black hole mass of at least several billion solar 
masses. A virial estimate of the black hole mass of 3 x IO'Mq was 
obtained by Willott et al. (2003) using the Mgn /12800 line-width 
and the scaling relations derived by McLure & Jarvis (2002). Us- 
ing the same scaling relations but also the virial estimator provided 
by Vestergaard (2002), Barth et al. (2003) consistently found a BH 
mass in the range (2 - 6) x 10''Mo from both the widths of the Civ 
/11549 and Mgn emission lines. Beside the factor 2.5-3 estimated 
uncertainty in the adopted scaling relations, this method assumes 
that the dynamic of the line-emitting gas is dominated by gravi- 
tational forces. Note, however, that for objects radiating close to 
the Eddington limit, the effect of radiation pressure cannot be ne- 
glected and may lead to a non negligible upward correction of the 
estimated BH masses (Marconi et al. 2008; Lamastra et al. 2010). 

2.2 Host galaxy properties 

High-resolution VLA observations of J 1148 have enabled to re- 
solve CO (3-2) line emission both spatially and in terms of velocity 
(Walter et al. 2004). The molecular gas is extended to a radius of 
~ 2.5 kpc from the center and a total mass of Mh, = 1.6x lO'^Mo 
has been derived, assuming a CO intensity-to-gas mass conversion 

' Schmidt, Schneider & Gunn (1995) find Mb = M1450 + 
2.5Q-Log(4400A/1450A) + 0.12 where the factor 0.12 comes from 
the zero point difl'erence between the AB and the Vega-based magnitude 
system for QSO-like spectra. 



factor appropriate for ultra luminous infrared galaxies (ULIRGs, 
Walter et al. 2004). Assuming that the gas is gravitationally bound 
and forms an inclined disk with inclination angle Walter et al. 
(2004) infer a dynamical mass of Miy^sin-i ~ 4.5 x 10'°Mg. The 
molecular disk inclination angle is the main source of uncertainty 
in the dynamical mass estimates (Wang et al. 2010). Walter et al 
(2004) correct for an inclination of 65°, on the basis of the broad 
CO FWHM detected, obtaining Mdy„ ~ 5.5 x IO'^Mq with 50% 
uncertainty. 

The derivation of the stellar bulge mass is extremely uncer- 
tain. Indeed this can be either estimated by extrapolating the local 
Mbh - M, relation (Marconi & Hunt 2003; Haring & Rix 2004) 
at high redshift, Mbh/0.0014 ~ 2.14 x IO'-Mq (if a BH mass of 
3 X 10' Mq is considered), or by subtracting from the observed dy- 
namical mass the H2 gas mass, M, ~ Mjyn - Mh, ~ 3.9 x 10'" MO. 
It is evident that these estimates are unconsistent, since they differ 
by almost two orders of magnitude. Such a discrepancy persists, 
though alleviated, even if we assume that the stellar bulge mass of 
the host galaxy of Jl 148 was 10 times the scale length of the molec- 
ular gas emission. Assuming a density profile of p oc r"^, the stellar 
mass would be M, ~ 3.9 x 10" Mq (Walter et al. 2004). 

Further constraints on the evolution of the QSO host galaxy 
come from observations of the gas phase metallicity and elemen- 
tal abundances in both the broad and naiTow line regions (BLRs, 
NLRs). BLRs metallicity tracers are emission lines ratios, such 
as Feii/Mgii (e.g. Barth et al. 2003), Nv/Civ (e.g. Pentericci et al. 
2002) and (Siiv+Oiv)/Civ (Nagao et al. 2006; Juarez et al. 2009), 
together with other metal lines like Cii (e.g. Maiolino et al. 2005), 
and Oi (Becker et al. 2006). The observed ratios are similar to those 
observed in low-z QSOs, indicating BLRs metallicities which are 
several times solar, ~ 7 Z© (Nagao et al. 2006, Juarez et al. 2009) 
and suggesting that the ISM of the QSO host was significantly en- 
riched by heavy elements. However, the BLRs are small nuclear 
regions, of less than a few pc, containing a total mass of ~ 10''Mo; 
thus, we consider such a high metallicity as an upper limit. Con- 
versely, NLRs have larger masses and a size which is roughly com- 
parable to the size of the host galaxy and represent a good tracer of 
the chemical properties on galactic scales, since their metal content 
reflects the past star formation history of the galaxy (Matsuoka et al. 
2009 and references therein). To date however, only measurements 
of the NLRs metallicity up to redshift 4 are available. Focusing on 
the Civ/Heii and Ciii/Civ flux ratios, Matsuoka et al. (2009) point 
out that even NLRs show no significant metallicity evolution with 
redshift, up to z ~ 4. For the highest redshift sample these authors 
find a NLR metallicity of about 1.32+[ j^Zq. We adopt this value as 
a lower limit for the ISM metallicity of Jl 148. 

2.3 Dust mass, FIR luminosity and SFR 

The dust present in the host galaxy is probably characterized by 
multi-temperature components: hot (~ 100 K), warm (30 - 80 K) 
and cold (10-20 K). The hot dust component, directly heated by 
the AGN activity, has been observed for several z ~ 6 QSOs with 
Spitzer (Charmandaris et al. 2004; Jiang et al. 2010) and repre- 
sents the least contribution to the total mass of dust. Warm dust, 
associated to starburst regions, dominates the emission in the rest- 
frame far-IR. For Jl 148, it has been observed in 4 bands: at 1.2 mm 
(Bertoldi et al. 2003), SSOfim, 450//m (Robson et al. 2004), and 
350//m (Beelen et al. 2006). The cold dust, finally, may represent 
the largest component. Its emission, however, would be hidden by 
the warmer component even if its mass were 1.5-3 times larger 
(Dunne et al. 2000). 
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In the rest-frame far-IR the emission is optically thin, T^iA) « 
1, and the dust mass can be estimated from the flux observed in a 
given band S yg as, 

Md„„ = ° ^ , (1) 

(l+z)K,(v)B(v,T,y 

where A:d(v) is the opacity coefficient per unit dust mass, B(v, T^) 
is the Planck function for a dust temperature Tj, and d]^ is the lu- 
minosity distance to the source. In the Rayleigh-Jeans part of the 
spectrum, dust radiates as a 'grey-body' with K^iv) = Ko(v/vo'f. In- 
tegrating over the spectrum one obtains the corresponding far-IR 
luminosity, 

Lfir = 4;rMdi,st J Ki(v) B(v, Ti) dv. 

Although in principle with four photometric points it could be pos- 
sible to simultaneously fit Mjust, yS, and Tj, at the redshift of Jl 148 
the observations do not fully sample the Rayleigh-Jeans part of the 
spectrum, preventing a determination of fj. Therefore, the estimated 
mass of dust and FIR luminosities for Jl 148 quoted in the literature 
have been derived under different assumptions for /3 and kq. The 
resulting dust masses, FIR luminosities and dust temperatures are 
listed in Table[T]and have been calculated through ax^ fit using all 
the observed points in the range 350 fjm - 1200 mm. Taking into 
account the intrinsic error on the estimated fluxes, the typical error 
associated to a single fit is 25%. As it can be inferred from the Ta- 
ble, we have considered, in addition to dust absorption coefficients 
already adopted in previous works (cases a-c), a fit to the Weingart- 
ner & Draine (2001) model for the SMC in the relevant frequency 
range (case d) as well as the dust optical properties expected for SN 
dust following the model by Bianchi & Schneider (2007, case e). 

The resulting dust masses and temperatures are very sensitive 
to the adopted dust absorption coefficient and spectral index with 
values in the range [1.86 - 4.78] X 10** Mo. In what follows we will 
assume an average dust mass of 3.4 x 10* with errorbars com- 
puted assuming the minimum and maximum values quoted in the 
Table. It is important to note that a mixture of SMC and SN ex- 
tinction curves (30% SMC and 70% SN) is found to best fit the 
extinction properties of a sample of QSOs with 3.9 < z < 6.4, 
including Jl 148 (Gallerani et al. 2010). 

From Table [T] we also infer an average FIR luminosity of 
LpiR = 2.2 X 10'^^ Lq with 15% uncertainty associated to the fit. As- 
suming that the dominant dust heating mechanism is radiation from 
young stars, it is possible to estimate the star formation rate (SFR) 
using the relationship between FIR luminosity and SFR derived by 
Kennicutt (19980, 

SFR _ LpiR 
IMoyr-i ~ 5.8 X 10' Lq' 

which implies a SFR ~ (3.8 ± 0.57) x W^MJyr (Bertoldi et al. 
2003). Such a high SFR is also supported by the first detection of 
the carbon [CII] line at 158 fim (Maiolino et al. 2005), where a SFR 
~ 3 X lO^Mo/yr is estimated combining the above conversion with 
the [CII] line - FIR luminosity ratio obtained from PDR models. 

It is clear that if the AGN contributes to dust heating, the cor- 
responding SFR would be much lower, consistent with the value 
that would be inferred from the observed gas mass by a simple ap- 
plication of the Schmidt-Kennicutt law. In fact, using the relation 

^ The conversion factor 5.8 X lO' has been derived assuming a mean lumi- 
nosity for a 10-100 Myr burst, stars with solar abundances and a Salpeter 
IMF (see the original Kennicutt 1998 paper for details). 



proposed by Daddi et al. (2010) to fit local ULIRGs and submil- 
limetre galaxies (SMGs) and QSOs, 

LogEsFR/[Moyr-'kpc--] = 1.42 x LogE„,J[Mopc--] - 2.93 

and assuming that the SFR is confined in the same central 2.5 
kpc where molecular gas emission has been detected, we find a 
SFR ~ 180Mo/yr consistent with previous findings by Dwek et al. 
(2007) and Li et al. (2007). 

Summarizing, the observed/inferred properties of J 1148 used 
as reference values in this work are: 

• A DM halo mass of Mi, = 10'^ Mq, which is assumed to host 
the z = 6.4 QSOJ1148. 

• A SMBH mass of Mbh ~ S^fjj x 10' M^, inferred from the 
Mgii ^2800 and dv ^1549 line-widths (Willot et al. 2003; Barth et 
al. 2003). 

• A dynamical and molecular gas masses of M^yn = (5.5 ± 
2.75) X 10'" Mo and Mh, = 1.6 x 10'° Mo, derived from obser- 
vations of the CO(3-2) line emission and assuming an inclination 
angle / = 65° (Walter et al. 2003). 

• A final steflar mass ranging from M, = (3.9 ± 2.75) x 10'" Mo 
(computed as the difference between the dynamical and molecular 
gas masses given above) up to M, ~ 10'^ Mq which is the one 
inferred if the local Mbh " M, relation were to hold even at redhift 
6.4. 

• An ISM metallicity of Z = I.32| Zo, estimated from obser- 
vations of NLRs in high redshift QSOs (Matsuoka et al. 2009). We 
refer to this value as a lower limit for the metallicity of Jl 148. 

• A total mass of dust of Mjus, = 3.4+] x lO** Mo, computed 
using eq. l[T) and the parameters given in Table[T] 

• An upper limit to the SFR, ~ (3.8±0.57)x 10' Mo/yr, obtained 
assuming that dust heating is entirely due to radiation from young 
stars. 



3 MODEL DESCRIPTION 

Here we describe our Monte Carlo semi-analytical model, GA- 
METE/QSOdust, which follows the formation and evolution of the 
quasar SMBH and of its host galaxy in the hierarchical scenario 
of structure formation. In this framework, a galaxy form through a 
series of merging episodes of lower mass fragments, called progen- 
itors. The code is divided into two main blocks: the first block runs 
backward in time and it is needed to reconstruct the merger tree 
history of the assumed lO'^Mo DM halo atz = 6.4 (see section|2]l; 
the second block runs forward in time and allows us to follow the 
evolution of the baryonic component and the growth of the nuclear 
black hole. In this way, it is possible to simulate several hierarchical 
merger histories of the host DM halo and then follow the build up 
of the quasar and the evolution of the host galaxy properties along 
these various formation paths. Each formation path will result in 
a peculiar redshift evolution for the host galaxy properties. How- 
ever, in order to draw some general conclusions, in what follows 
we will show the results obtained averaging over all the simulated 
hierarchical merger histories. 

In the following sections we will briefly describe the Monte 
Carlo algorithm used to reconstruct the hierarchical merger tree 
of the lO''' Mo DM halo hosting the QSO at redshift z = 6.4, 
(see SSF07 for a detailed description of the algorithm) and the 
semi-analytical model which follows the evolution of gas, stars 
and metals along the hierarchical build-up of the galaxy (see Sal- 
vadori, Ferrara & Schneider 2008, herafter SFS08, for more de- 
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Table 1. Dust temperature, dust mass, and FIR luminosity resulting from the^^-fit of the four observed photometric points of Jl 148 (see text). Each lines refer 
to the values inferred from the fit taking the dust absorption coefficient as specified by the corresponding atq and j8. These values have been chosen following 
quoted values in the literature and assuming specific dust grain models. In particular: a = Bertoldi et al. (2003); b = Robson et al. (2004); c = Beelen et al. 
(2006); d = fit to the Weingartner & Draine (2001) model for the SMC in the spectral range [40-200];jm; e = SN dust model by Bianchi & Schneider (2007). 



Ref. 


Ko [cm^/gr] 






Td [K] 


Mdust [Mq] 




a 


7.5 


230 


1.5 


58 


3.16 X 10** 


2.32 X 10'3 


b 


30 


125 


2.0 


49 


2.91 X 


2.09 X 10'3 


c 


0.4 


1200 


1.6 


56 


4.29 X lO** 


2.27 X lO''' 


d 


34.7 


100 


2.2 


47 


4.78 X 10** 


2.02 X 10'3 


e 


40 


100 


1.4 


60 


1.86 X 10** 


2.38 X 10'3 



tails). Then, we will introduce in details the new features of the 
GAMETE/QSOdust version implemented to follow the evolution 
of the central SMBH and of the dust mass. 



3.1 The Hierarchical Merger History 

The possible hierarchical merger histories of a M;, = 10'^ Mg 
DM halo at redshift z = 6.4 are reconstructed using a binary 
Monte Carlo algorithm with mass accretion (e.g. Cole et al. 2000; 
Volonteri, Haardt & Madau 2003), based on the Extended Press- 
Schechter (EPS) theory (e.g. Lacey & Cole, 1993). 

At each increasing redshift, a DM halo can either loose part 
of its mass or loose mass and fragment into two progenitors with 
random masses in the range Mres < M < Mi, 12, where M^^s is a 
threshold mass representing the model resolution limit: halos with 
M < Mres, are not considered as progenitors but account for the 
surrounding medium in which other halos are embedded and from 
which they accrete mass and virialize. 

The resolution mass, M^s, along with the redshift interval 
chosen to sample the merger history, represent the free parame- 
ters of the model and are constrained in order to ( i) prevent mul- 
tiple fragmentations, (ii) resolve all the progenitors with masses 
above a given physical threshold, (Hi) reproduce the EPS predic- 
tions, (iv) control the computational cost. Given the above require- 
ments, we assume that Mfes(z) = lOM(Ty,y = lO^'K, z), where 
M(rv,r = 10'*K, z) is the mass of a DM halo with a virial equilib- 
rium temperature of lO^'K at redshift z- The merger tree is computed 
using 5000 redshift intervals logaritmically spaced in the expansion 
factor between z = 6.4 and z = 37. In Fig. [T] we compare the mass 
function of the progenitor halos averaged over 50 possible hierar- 
chical merger histories with the EPS predictions. We can see that 
the agreement is extremely good at all redshifts (z= 15, 13, 10, 7). 

We then use the produced merger trees as an input to recon- 
struct the formation and evolution of the quasar Jl 148 and its host 
galaxy through cosmic time. 



3.2 Evolution along a merger tree 

In what follows we present a brief summary of the basic prescrip- 
tions proposed by SFS08 which are unchanged in the present ver- 
sion of the code. Note that in order to significantly reduce the high 
computational cost caused by the high number of redshift intervals, 
we are now considering as star-forming objects all DM halos with 
masses equal or greater than M^fiz) = 10M|.es(z) (see Fig.O. 

In any star-forming halo, the star formation rate is taken to be 
proportional to the mass of cold gas. 
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Figure 1. Number of progenitors of a QSO host halo of mass 10 M© at 
z = 6.4 as a function of mass at four different redshifts: z = 15, 13, 10 and 7 
(from top left to bottom right). In each panel, histograms represent averages 
over 50 realizations of the merger tree and the errorbars indicate Poissonian 
error on the counts in each mass bin; solid lines show the predictions of EPS 
theory and the vertical lines indicate the values of the resolution mass at the 
corresponding redshift. 



whose gradual accretion is regulated by an infall rate, dM,„f/dt 
(see SFS08 for the complete expression). The parameter f,(z) is 
the redshift-dependent global efficiency of star formation per unit 
time. We will discuss this efficiency in detail in the next section. 

In this model, Population III (Pop III) stars with an average 
mass of 2OOM0 are assumed to form if the gas metallicity is lower 
than a critical value = 10"''Zo. Otherwise, if the gas metallicity 
exceeds the critical limit, low-mass (Pop II/I) stars are assumed to 
form according to a Larson initial mass function (IMF), which fol- 
lows a Salpeter-like power law at the upper end but flattens below 
a characteristic stellar mass: 



0(m) oc m 



(3) 



SFR = /.(z)Mga„ 



(2) 



where a = 1.35, m^, = 0.35Mg and we normalize the integral of 
m(j)(m) in the mass range (0.1 - 100) Mq to unity. The evolution 
of very massive stars is rapid (they reach the end of their main 
sequence phase in 3-5 Myr) and ends in a violent explosion as pair 
instability supernovae (PISN) which leave no remnants behind and 
significantly contribute to the metal enrichment of the interstellar 
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medium at the lowest metallicities. Thus we can assume that Pop 
III stars die instantaneously, while less massive stars are assumed to 
evolve according to their appropriate stellar lifetime. In this work, 
the lifetimes of stars of different mass and metallicity are computed 
according to the parametric form proposed by Raiteri et al. ( 1996). 
As we will show in section [6.3.21 our results are not sensitive to 
the adopted values for the mass of Pop III stars and of the critical 
metallicity. 

Finally, the galactic outflows driven by SN explosions, 
dMej^sti/dt, are assumed to be proportional to the SN explosion rate 
(Rsn) divided for the escape velocity (v^) squared: 



dm, 



■cj.SN 



dt 



(4) 



where 6,^. is the SN-wind efficiency and £sn is the average SN ex- 
plosion energy (see SFS08 for further details). 
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3.3 New features 

The new features implemented for the purposes of the present work 
are: (i) the growth of the BH (via both gas accretion and mergers 
with companion BHs), ( ii) the effect of the AGN feedback on the 
host galaxy, (Hi) the merger driven starbursts, i.e. bursts of star for- 
mation episodes enhanced by halos binary merging, and (iv) the 
dust formation and evolution into the ISM. We will describe them 
in turn. 



Figure 2. The threshold masses adopted in the code as a function of red- 
shift. M(ltf'K, z) is the mass of a halo with a virial temperature Tyif = lO'* 
K, M,-es(z) = 10M(Tyij, z) is the minimum halo mass in the merger tree sim- 
ulations and Msf (z) is the minimum halo mass to form stai's, assumed to be 
= 100M(10*K,z). These masses are compared with the 4-0" mass density 
fluctuation M4_£r(z), green solid line. When newly virialized halos, in the 
merger tree, are less massive than M4_o-(z) they do not host a seed black 
hole. The redshift at which halos are no more populated with seed BHs de- 
pends on the assumed threshold as indicated by the vertical lines. For the 
adopted minimum halo mass to form stars, M^f (z), we populate halos down 
to redshift of about 8. 



3.3.1 BH growth 

It is usually assumed that SMBHs are assembled by mergers with 
other BHs and/or by accretion of gas from the surrounding medium. 
Thus, the growth of a BH is determined by the seed BHs properties 
(e.g. mass, redshift distribution and abundance), the efficiency of 
the accretion rate, the number of BH mergers and the dynamics 
of the coalesced BHs. Theoretical studies suggest that a seed BH 
with mass in the range [10^ - 10*]A/o can form either by the rapid 
collapse of Pop III stars (Heger & Woosley 2002) or by the direct 
collapse of massive hot (T > lO^'K) and dense gas clouds induced 
by gravitational instabilities (e.g. Bromm & Loeb 2003; Begelman, 
Volonteri & Rees 2006; Volonteri & Natarajan 2009). 

In this work we assign a BH seed, Msecd = IO^'/i^'Mq to all 
the progenitor halos corresponding to density fluctuations higher 
than 4 - cr {M > M^_^) that exceed the threshold mass required 
to form stars {M > Mjf) for the first time in the merger tree. The 
evolution of these critical masses as a function of redshift is shown 
in Fig|2l From this plot we can infer that the probability, P, that a 
halo with M > M^i will host a seed BH depends on the redshift, 
other than on the specific formation history, and that P — » for 
z < 8. Note that even assuming P = 1 for all the star forming 
progenitors, the mass density parameter in seed BHs at redshift z = 
6.4 would be Q., ~ 10"", several orders of magnitude less than the 
SMBH density parameter estimated from nuclear BHs of nearby 
galaxies, nsMBH(z = 0) = 2.9 x lO"*^ (Merloni & Heinz 2008). 
Indeed, seed BHs grow via binary mergers and gas accretion. The 
choice of the BH seed will be discussed in section [Ol Following 
Tanaka & Haiman (2009), we assume that if the progenitor halos 
mass ratio, = M^/Mp (where is the mass of the secondary 
smaller halo and Mp is the mass of the primary larger halo), is yu > 
1/20 the BHs coalesce elficiently over the same timescale of their 



host halos, adding their masses linearljU Conversely, if the mass 
ratio of the two progenitors is smaller than 1:20, only the largest of 
two original BHs survives as a nuclear BH, the less massive one is 
assumed to end up as a satellite and its evolution along the merger 
tree is no longer followed. 

It is known that the BH produced by a coalescence receive a 
gravitational recoil (kick velocities > 100 km/s) due to the net lin- 
ear momentum accumulated by the asymmetric gravitational wave 
emission (Campanelli et al. 2006; Baker et al. 2006). As a con- 
sequence of this effect, the recoiling BH could be kicked out of its 
host halo or displaced in under-dense regions, thus reducing the gas 
accretion rate. Tanaka & Haiman (2009) have added in their Monte 
Carlo code an explicit calculation of the orbits of kicked BHs, and 
self-consistently include their corresponding time dependent accre- 
tion rate. We do not attempt to do a similar analysis here and we 
neglect this effect; note, however, that in our model we expect the 
number of BHs ejected out of their host halo to be modest since the 
escape velocities of halos with masses > Msf(z) are > 1 10 km/s. 

In GAMETE/QSOdust, BHs are allowed to grow in mass by 
accreting gas from the surrounding medium. This process is trig- 
gered by gas-rich halo mergers and it is regulated both by star 
formation, which consumes gas by converting it into stars, and by 
AGN feedback, which reduces the amount of material available for 
accretion. In our formulation, the central BH accretes gas at a rate 
given by M^cci = ™'"(^Edci> ^bhl)^ where Mem is the Eddington 
rate, 

4nGMBHinp 

Mem = (5) 

ercrjc 

' This condition is derived requiring that the halo merger timescale is 
shorter than the Hubble time (see Tanaka & Haiman 2009 for more details). 
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and Mbhl is the Bondi-Hoyle-Lyttleton (BHL) accretion rate, 
Mrhi. - ^ ■ 



(6) 



In the above equations, G is the gravitational constant, lUp is the 
proton mass ctt is the Thomson cross section, c is the light speed, 
Cs is the sound speed, Pgas('"A) is the gas density evaluated at the 
radius of gravitational influence of the BH (or Bondi accretion ra- 
dius), = 2GMbh/Cs, and or is a dimensionless parameter (see 
below); in eq. |5] the parameter is the radiative efficiency, which 
determines the conversion efficiency of mass accretion into energy 
released as radiated luminosity, and it is fixed to 0.1 (Shakura & 
Sunyaev, 1973). 

We adopt a gas density distribution described by a singular 
isothermal sphere (SIS) profile with a flat core: 



p(R) 



1 + (R/R,„,,)^ 



(7) 



where R^orc = 0.012i?vir is the size of the flat core, i?vii is the virial 
radius and p„omi is the normalization constant; the profile is nor- 
malized such that, at each time step, the total gas mass, M^as, is 
enclosed within the virial radius. The factor a which appears in 
eq. [6] does not appear in the original analysis of Bondi & Hoyle 
(1944) and Hoyle & Lyttleton (1939), but has been introduced by 
Springel, Di Matteo & Hernquist (2005) as a numerical correction 
factor, to compensate for the limitations of the numerical simula- 
tions. In fact, present simulations lack both the resolution and the 
physics to model the multiphase interstellar medium, and tend to 
strongly underestimate the Bondi-Hoyle accretion rate; values of 
100 < Q- < 300 are commonly adopted in the literature (Springel 
et al. 2005; Di Matteo et al. 2005, 2008; Sijacki et al. 2007; see 
however Booth & Schaye 2009 who parametrize o- as a function 
of density). As it will be discussed in the following sections, we 
also assume a constant a and fix its value in order to reproduce the 
observed black hole mass of Jl 148. 



3.3.2 AGN feedback 

It is natural to expect that a quasar shining at, or close to, its lim- 
iting Eddington luminosity can generate a powerful galactic wind 
and eventually terminate the accretion process that feeds it. If the 
energy in the outflow liberates as much energy as the binding en- 
ergy of the gas in a dynamical time, the feedback from the ac- 
creting BH can be self-regulated (Silk & Rees 1998; Wyithe & 
Loeb 2003). Hydrodynamical simulations suggest that BH accret- 
ing at high rates release enough energy to drive galactic scale out- 
flows of gas (Springel et al. 2005; Di Matteo et al. 2005; Ciotti, 
Ostriker & Proga 2009, 2010). The mechanism by which the en- 
ergy released by the accreting BH is coupled to the surrounding 
medium is currently under debate (see Ciotti et al. 2010 and refer- 
ences therein). The energy transfer rate is commonly parametrized 
as (e.g. Springel et al. 2005), 



■ fw.AGNfrA^ac 



(8) 



where the parameter 6„_agn is the coupling efficiency and it is usu- 
ally assumed to be a free parameter of the simulation, independent 
of the environment and gas properties. 

Note that recently, Sijacki et al. (2007) have dilferentiated be- 
tween a 'quasar mode' and a 'radio mode' feedback, depending 
on the accretion rate. In particular, BH-powered winds are related 
to the luminous quasar activity and are efficient for high accretion 
rates. They are expected to be more powerful than starburst- driven 



winds and should be able to sweep away a significant amount of gas 
from the deep potential well of their massive host galaxies. Once 
the obscuring material is swept away from the center of the system, 
the quasar enters its visible phase in which it outshines the entire 
galaxy (quasar mode). In the low-accretion rate regime, which is 
subdominant in terms of black hole mass growth, the feedback is 
characterized by AGN-driven bubbles; this regime can be identied 
with radio galaxies in clusters (radio mode). The BH accretion rate 
which signs the transition between these two modes is usually as- 
sumed to be (Miccr/MEdd) = 10"^ (Sijacki et al. 2007); thus, 
since we are not interested here in the radio-mode regime, we only 
activate the BH feedback when the accretion rate is super-critical. 

When the above condition is verified, we estimate the amount 
of gas mass ejected out of the host halo per unit time as. 



liMej.AGN , / C \2 
; = 2ew,AGNei — M^, 

at ^ V, ' 



(9) 



where Mej,AGN is the mass of ejected gas and Ve = (f^) is the 
escape velocity of the galaxy. 



3.3.3 Bursting star formation mode 

The star fomation history (SFH) of the QSO host at redshift z = 6.4 
is the result of the SFHs of all its progenitor halos. There are ob- 
servational and theoretical indications that support the idea that 
galaxy mergers have a dramatic effect on their morphological evo- 
lution and star fomation history (see e.g. Woods & Geller 2007; 
Ellison et al. 2010; Mihos & Hernquist 1996; Cox et al. 2008 
and references therein). In particular, enhanced central star forma- 
tion in interacting galaxies result from gaseous inflows that occur 
when the gas loses angular momentum through gravitational tidal 
torques produced primarily by the non-axisymmetric structure in- 
duced by the companion galaxies. Indeed, there is mounting evi- 
dence that near equal mass mergers (with mass ratios greater that 
1:3) triggers the most vigorous star-forming galaxies in the local 
universe, the ULIRGs (see e.g. Genzel et al. 2001; Tacconi et al. 
2002; Dasyra et al. 2006). Similarly, there is growing observational 
evidence for interaction-driven nuclear activity (see the recent re- 
view by Shankar 2009). 

Based on a set of numerical simulations. Cox et al. (2008) 
have studied the effects of the galaxy mass ratio onto merger-driven 
starbursts. They find that while mergers between nearly equal mass 
galaxies produce the most intense bursts of star formation (major 
mergers), minor mergers produce relatively little enhancement with 
respect to the quiescent mode of star formation. These simulations 
have been calibrated to reproduce the observed properties of lo- 
cal disk galaxies and may not be immediately extended to higher 
redshifts. The resulting fitting formula have been implemented in 
semi-analytic models of galaxy evolution which successfully re- 
produce a number of observations at z ~ (Somerville et al. 2008). 

Here we are interested in exploring the effects of different 
SFHs on the chemical properties of QSOs host galaxies. There- 
fore we investigate both quiescent and bursting models: in the first 
class of models, star formation occurs only in a quiescent mode 
which is independent on the galaxy mergers mass ratio; conversely, 
in the second class of models an additional contribution to the stel- 
lar mass comes from galaxy major mergers. In order to explore a 
wide parameter space, we parametrize the global redshift depen- 
dent efficiency of star formation per unit time that appears in eq. |2] 
as. 



f,{z) = (fquics + eburst)/fdyn(z) 



(10) 
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where ?(jyn(z) = R^iJvc is the dynamical time, equie, is the constant 
star formation efficiency during the quiescent phase of the evolu- 
tion and 6buist is the starburst efficiency, expressed as a normalized 
Gaussian distribution of the mass ratio. 



1 



^buist 



^/2nt 



(11) 



urst 



The starburst efficiency peaks at the critical mass ratio pciit = 1 
and the dispersion parameter controls the intensity of the burst for 
a given mass ratio. It is important to stress that since the merger 
rate and the mass ratio of merging halos depend on the hierarchical 
merger history, variations of the parameter tr allow to explore 
a wide range of SFHs, with lower values of tTburst coiTesponding 
to fewer intense bursts over the history of the galaxy; as iTburst de- 
creases, the SFHs will be characterized by more frequent and less 
intense bursts, progressively approaching the quiescent SF mode. 
For our test cases, cTbmst = 0.05,0.1,0.25 corresponding to maxi- 
mum mergers efficiencies (if /j = /Jcdt) of ei,„^^t - 8,4, 1.6, respec- 
tively. 



3.3.4 The evolution of dust 

We compute the evolution of dust under the following assumptions: 

(i) both SNe and AGB stars contribute to dust production on dif- 
ferent timescales according to their stellar progenitors lifetimes; 

(ii) there is no dust infall or outflow to/from the ISM, i.e. the in- 
falling gas is dust free: only dust grain metals which are returned to 
the gas-phase ISM through destruction by interstellar shocks (i.e. 
dust which has been shocked by SNe) can be ejected out of the 
galaxy; (Hi) a fraction f„,c of the dust mass is enclosed in dense 
MCs, M^it) = fmcMj(t); this component can grow through mass 
accretion of metals onto pre-existing grains in a time scale t^^c', f'^J 
the remaining dust component, M''j^^ {t) = (1 - fmc)Mj(t), is dif- 
fused in the ISM and can be destroyed by interstellar shocks in a 
time scale t^. 

The destruction and accretion time scales depend on the prop- 
erties of the ISM and on the cycling times between the cloud (cold 
and dense) and intercloud (waim/hot and diffuse) phases (e.g. Tie- 
lens et al. 1998, 2005). Thus, a multiphase description of the ISM 
would be required to properly follow the evolution of the two dust 
components, M""^ and Mj'-^-^. Here, we make the simplifing assump- 
tion that a fixed fraction f„,^ ~ 1/2 (see e.g. Dwek 1998) of the total 
dust particles are located into MCs and thus can serve as growth 
centers for the accretion of icy mantels (e.g. Draine 1990). 

Following Hirashita (1999), we model the dust accretion 
timescale T„,.r as: 



(12) 



^row[l - fdepit)] 



where, T„,r = 3 x lO' yr is the typical lifetime of a MC (e.g. Tielens 
1998), X„ij is the fraction of cold gas into the ISM, which is set 
to unity in our one-phase scenario. Finally, fdepit) is the depletion 
factor, i.e. the ratio between the mass of dust and the total mass in 
metals Mj(t)/Mz(t) which is computed at each time in each halo 
of the simulation. Only metals which are not already condensed 
into dust grains (i.e. gas-phase metals) are available for accretion. 
Typical values of this timescale are in the range (10^ - 10**) yr (see 
Fig.|9). 

The destruction timescale, r^, i.e. the lifetime of dust grains 
against destruction by thermal sputtering in high- velocity (v > 150 
km s~') SN shocks is given by (e.g. Dwek & Scalo 1980; McKee 
1989) : 



Td 



Mism(0 



(13) 



where is the dust destruction efficiency, Ms„ep, is the effective 
ISM mass that is completely cleared of dust by a single SN rem- 
nant and R'g^ = /sn^^sn is the effective SN rate for grains destruc- 
tion, which takes into account that not all the SNe which inter- 
act with the ISM are efficient in destroying dus0 (see e.g. McKee 
1989; Tielens 2005); for PISNe fsN = I and for core collapse SNe 
.fsN ~ 0.15 (i.e. 15% of all SNe in a stellar generation are efficient). 
Following McKee (1989), we compute the mass of the interstel- 
lar matter shocked to a velocity v,;, in the Sedov-Taylor stage, in a 
homogeneous ISM, as: 



6800Mp 



(£sN/10^'erg) 
'(v,,,/100kms-')2 



(14) 



where v,;, ~ 200 km s"' is the (minimum) non-radiative shock ve- 
locity and EsN is the average supernova explosion energy, assumed 
to be 2.7 X 10^^ erg for PISNe and 1.2 x lO""' erg for core col- 
lapse SNe. We adopt the dust destruction efficiencies from Nozawa 
et al. (2006) for PISNe and from Jones et al. (1996) for SNe. As- 
suming a reference value for the ISM density of /iism = 1cm"' 
we get frf ~ 0.6(0.48) for PISNe (SNe). With these parameters, 
we find that for conditions which apply to the Milky Way Galaxy 
{M,sM ~ 5 X 10' Mo, R'^f^ ~ 7.5 x 10"' yr"', McKee 1989), typical 
grain lifetimes are ~ 0.6 Gyr, in agreement with Jones et al. (1996). 



4 THE CHEMICAL EVOLUTION MODEL WITH DUST 

In this section we describe the chemical evolution model imple- 
mented in the code GAMETE/QSOdust. 

We will indicate with Mjsm the total mass in the interstellar 
medium (Mg^s + Mj) and with Mz the total mass of metals (dif- 
fused in the gas phase and condensed into dust grains). At any given 
time, the mass of gas phase metals is calculated as the difference 
between Mz and Mj (with proper stoichiometric coefficients when 
we are interested in a particular gas phase element). The variations 
of the total mass of gas, stars, metals and dust is followed in each 
progenitor halo of the hierarchical tree solving the following set of 
differential equations: 



dM,(t) 
dt 



: SFR(f) ■ 



dR(t) 
dt 



dt 



_SFR(0 + ^ + '^^-f^'^ '^^''^^^ 



-d - f.) 



dt 

dM„,„{t) 
dt 



dt 



dt 



dMzit) dYzU) dMi„f(t) 

- -Zism(OSFR(0 + + Z„,(f)- ' 



dt 



dM,j(t) 
dt 



dt 



-ZisM(f)(l-e,) 



dt 

dMa,„{t) 
dt 



(15) 



(16) 



(17) 



The most massive stars exploding in an association will be very effective 
in destroying dust while sweeping up matter into a shell. Subsequent su- 
pemovae will interact predominantly with the gas inside the shell and the 
shell itself and, therefore, do not contribute much to dust destruction. The 
expansion of the supershell is too slow to destroy much dust (e.g. Tielens 
2005). 
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dt 



-Zrf(OSFR(?) + 



dt 



.dM„. 



dt 



(18) 



where M, and are the total mass of stars and dust, Zism(0 = 
Mz(t)/Misu(t) is the ISM metallicity, Z,,(/) = Mj{t)/MisM(t) is the 
total dust abundance in the ISM and Zyj, is the metallicity of the 
infalling gas, namely of the hot gas at the virialization epoch. Its 
value corresponds to the metallicity of the external medium, that is 
the ratio between the total mass of metals ejected by star forming 
haloes and the mass of the diffused gas, not enclosed in collapsed 
objects (see SFS08 for details). 

The new features of the code are summarized by: (i) the evo- 
lution of dust described by eq. [18] (ii) the term (1 - er)dM„ccr/dt 
which accounts for the fraction, (I - er), of the gas (metals and 
dust), which fuels BH growth via accretion; (Hi) the gas ejection 
rate dMej{t)ldt = dM^i^sNiO/dt + dMfjAGN(t)/dt which includes 
both SNe and AGN feedback. 

Finally, the terms dR/dt and dYz/dt and d¥j/dt are the rates 
at which gas, heavy elements and dust are returned to the ISM. 
The equations and grids of yields used to compute the gas and 
metal return fraction are the same as in SFS08 and V09: van den 
Hoeck & Groenewegen (1997) for AGB stars with initial metal- 
licities Z = (5 x 10"^ 0.2, l)Zo and masses (1 - 8)Mo; Woosley & 
Weaver (1995) for SNe with initial masses (12-40)Mo at metallici- 
ties Z = (0, 10-\ 10-^ 10^', I) Zo; Heger & Woosley (2002) for the 
PISNe assuming a reference progenitor mass of 200 Mq. The dust 
injection rate dYj(t)/dt is computed as in eq. (12) of V09, including 
contributions from AGB stars, SNe and PISNe. We use the grids of 
Zhukovska, Gail & Trieloff (2008) for AGB stars with O.SMa < 
m < 8M0 at different metallicities; for SNe with progenitor mass 
of (12-40)McE|andmetallicitiesZ = (0, 10-^ 10-^ 10^', l)Zo, we 
take the model by Bianchi & Schneider (2007) where grain con- 
densations starts from seed clusters made of A/' > 2 monomers with 
a sticking coefficient of o- = 1 . We also assume that SNe explode in 
an ISM with a density of pisM = lO^^^g cm"^, so that only 7% of 
the initial dust formed in the ejecta survives the passage of the re- 
verse shock. The resulting dust masses (10"'' - 10"' Mq depending 
on the progenitor star initial mass and metallicity) are in agreement 
with the values inferred from observations of young SN remnants 
(e.g. Hines et al. 2004; Rho et al. 2008; Sandstrom et al. 2008; Rho 
et al. 2009). Finally, for PISNe we adopt the dust yields computed 
by Schneider, Ferrara & Salvaterra (2004) for PISNe. As in V09 
we extrapolate from the dust mass produced by the largest avail- 
able AGB model (7 Mg) and by the smallest SN progenitor (12 
Mq) to account for metal and dust production in the mass range 
(8 - 11) Mg for which the dust yields are not yet available in the 
literature. The grids of Zhukovska et al. (2008) are not computed 
at metallicity Z < 5 x IO"^Zg; thus we have assumed that the dust 
mass produced by AGB at lower metallicities is the same as that at 
Z = 5 X lO-^zJl 



^ Dust and metal yields from Pop II/I stars with masses in the range 41 - 
lOOMg are assumed to be zero. Indeed, it has been shown (e.g. Dwek et al. 
1998; Zhukovska et al. 2008 and reference therein) that such massive stars 
provide only a minor contribution to the total dust in the ISM with respect 
to that of AGB stai's and core-collapse SNe. 

* Note that this approximation may hold for carbon dust but silicate dust 
production is likely to be significantly lower at low metallicity (see V09). 



5 MODEL PARAMETERS 

Our main aim is to investigate the effect of different SFHs in mod- 
eling the chemical properties of the host galaxy of QSO Jl 148. The 
free parameters which determine the shape of the SFH are: 

• the star formation efficiency, /,(z), which is determined by the 
choice of 6quies and (Xbuist (see eq.llOt; 

• the efficiency of the AGN-driven wind, e„ .AG;v (see eq.[9]l; 

• the efficiency of gas accretion onto the BHs, a (see eq.[6ll. 

The parameters which define the models considered in this 
work are summarized in Table |2] where letters and numbers given 
in the first column identify each model. 

These have been selected in order to reproduce the BH and 
total gas mass observed at z = 6.4. Given the uncertainties in the 
derivation of the stellar mass discussed in section [2^ we have not 
constrained evolutionary models in terms of the predicted final stel- 
lar mass. Indeed, by varying the global star formation efficiency, 
/.(z), we explore the dependence of the final mass of metals and 
dust on the stellar mass and SFH, considering models which pre- 
dicts a stellar mass compatible with the value inferred by the obser- 
vations (3.9 X lO'^Mg, Walter et al. 2004) as well as models with 
higher stellar masses, closer to the values implied by the present- 
day Mbh-A/. relation (~ lO'^Mg, see Fig.O. The rate at which gas 
is converted into stars in quiescent (Q) models is regulated by the 
choice of the efficiency eq„ies given in the second column of the Ta- 
ble (for these models fburst = in eq.lIOt. In bursted (B) models the 
frequency and amplitude of starbursts are regulated by the parame- 
ter cr buisi. The dependence on the choice of this latter parameter is 
investigated in Appendix A. 

It has been shown that assuming an AGN wind efficiency 
fw.AGN = 5% many observed properties of elliptical galaxies as well 
as the observed Mbh -o" correlation can be reproduced (see the dis- 
cussion in Hopkins et al. 2006). Here, we follow a similar approach: 
in all the models presented here, we choose the parameters a and 
SwAGN in order to reproduce the observed BH mass and to mimic 
a powerful AGN-driven wind which is able to sweep away enough 
gas from the host galaxy to reproduce the observed gas content and 
make the quasar optically visible, in its active phase, at z = 6.4. 
We set the parameter a in the range 180 - 350, to allow enough 
gas accretion to fuel the observed SMBH and the wind efficiency 
fvi.AGN, is chosen to regulate simultaneously the BH accretion, star 
formation and gas ejection (to produce a global gas mass in agree- 
ment with the observational lower limit of Walter et al. 2004, see 
section [2^ . 

Strictly speaking, the parameters listed in Table 2 are not in- 
dependent. With higher star formation efficiencies, /,, less gas will 
be available for BH growth, thus requiring a higher BH accretion 
efficiency, a, in order to reproduce the observed BH mass. On the 
other hand, the higher is BH accretion rate the stronger is the AGN 
wind (and the earlier it affects the global gas content), thus we need 
to set properly the coupling efficiency, e„ .AGN> to drive a gas outflow 
which expels the required amount of gas, but avoids a too early de- 
cline of the SFR and BH accretion. Therefore, as it has emerged 
from the extensive exploration of the parameter space out of which 
the present models have been selected, the set of parameters which 
characterize each model is not unique. Given this degeneracy, we 
keep the value e,i_AGN = 5 x 10"'' fixed as we vary a, in order to 
easily compare the AGN feedback effects in different models and 
control the combined effect of BH accretion and AGN-driven gas 
outffow, in shaping the SFH. Unfortunately, the AGN-driven out- 
flows are very difficult to observe if the quasar is in its active phase. 



10 R. Valiante, R. Schneider, S. Salvadori, S. Bianchi 



model 




C"buist 


^biirst,max 


a 


£vr,AGN 


Qi 


0.045 






190 


5 X 10-3 


Q2 


0.1 






180 


5 X 10-3 


Q3 


1.0 






350 


5 X 10-3 


Bl 


6.7 X 10"' 


0.05 


0.53 


180 


5 X 10-3 


B2 


0.02 


0.05 


1.6 


200 


5 X 10-3 


B3 


0.1 


0.05 


8.0 


200 


5 X 10-3 



Table 2. Model parameters. The models presented in this work are divided 
in two classes on the basis of the simulated SFH: quiescent only (Q) and 
bursting (B) star formation models. The numbers in the first column identify 
the efficiencies of star formation in quiescent phases. The second column 
provide the value of o"burst which determines the number, redshift distribu- 
tion and intensity of the starbursts (B-models only). The maximum starburst 
efficiency, eburst.max is given in the third column (B-models only). Finally, 
the efficiency of gas accretion into the central BH, a, and the efficiency of 
AGN-driven wind, ewAGN, are shown in the last two columns. 

As in the case of Jl 148, optically bright quasars outshine their host 
galaxy and ionize the cold gas which is usually a tracer of the out- 
flows. Thus, there are no observational constraints on the feedback 
process in such peculiar high redshift objects. 



6 RESULTS 

In this section we present the results of the semi-analytical hierar- 
chical model GAMETE/QSOdust. We have performed 50 random 
realizations of the quasar hierarchical assembly (as explained in 
section 3.1). Each of these merger histories allows us to predict a 
specific evolutionary path for the host galaxy chemical properties. 
However, to make the results and conclusions of this work inde- 
pendent of the hierarchical formation history, we will present the 
predicted evolution of the BH and galaxy properties averaging over 
all the 50 merger tree realizations. 



6.1 BH evolution 

As introduced in section |Z2] we investigate the efi'ects of the as- 
sumed SFH considering either models which predict stellar masses 
in agreement with the J 1148 position on the Mbh - M, plane and 
with what would be implied if the system were to follow the ob- 
served local relation, taking into account the observational disper- 
sion. Such models are shown in Fig. [3] which depicts the evolution 
onto the Mbh - ^* plane predicted by quiescent models Ql, Q2, 
Q3 (upper panel) and bursted models Bl, B2, B3 (lower panel). 
By construction, all models reproduce the observed J1148 SMBH 
mass but predict stellar masses increasing with /, : low-/, models 
(Ql and Bl) reproduce the M, inferred from observations; high-/, 
models (Q3 and B3) predict a final Mbh - M, correlation which is 
in agreement with the local observations (squared data points). We 
note that, despite the different final stellar masses, all the models 
tend to approach the final value from the bottom. In other words we 
find that in most of the models the stellar bulge appears to evolve 
faster than the central BH (M. > 10'° Mq always for Mbh > 10^ 
Mo). This differs from the recent study by Lamastra et al. (2010) 
who predict an evolution with redshift of the Mbh - ^» > with higher- 
z galaxies forming their nuclear BH faster than their stellar mass. 
This may be due to the different physical description of BH accre- 




IC IQO 109 101° 10'2 10'= 



Mst„ (Mo) 

Figure 3. The evolution of the BH mass as a function of the stellar mass pre- 
dicted by quiescent models (upper panel) and bursted models (lower panel). 
Open squared data points with error bars are from Marconi & Hunt (2003, 
Table 1 in their work) and represent the Mbh - M, relation observed in the 
local Universe, with the empirical fit Mbh/M, ~ 0.002 (dotted line) pro- 
vided by the same authors. The filled circle shows the Jl 148 BH (Willot et 
al. 2003) and stellar masses (Walter et al. 2004) inferred from observations. 
For all models, solid lines are averages over 50 merger tree reafizations and 
shaded areas represent 1 cr dispersion. 

tion and feedback which, in the Lamastra et al. (2010) study, tend 
to enhance the accretion efficiency at high redshift. 

For each model and at each redshift, we compute the total 
BH mass and accreted gas mass as the sum of all nuclear BHs and 
accretion rates hosted in each progenitor halcQ. These are shown 
in Fig.|4] Due to the strong correlation among the parameters dis- 
cussed above, the BH accretion rate (upper panels) and BH mass 
(lower panels) evolution have a negligible dependence on the par- 
ticular model. At redshift z = 6.4, we predict that J 1148 has a 
BH mass of ~ 3 x 10' Mq and accretes up to ~ (20-30) MJyv of 
gas; in units of the Eddington rate, the predicted accretion rates are 
MflH/MEdd ~ (0.3 - 0.45) for all the models investigated. The thin 
line in the BH mass evolution (bottom panels) represents the final 
BH mass that would be produced only through BH mergers. It is 
clear from the figures that even if BH mergers are the main drivers 
for the early growth of the BH mass at z > II, their contribution to 
the final BH mass at z = 6.4 is less than 1%. 

At redshift z > 1 1 the dependence on the particular hierarchi- 
cal formation history is stronger than at lower redshifts, as indicated 
by the larger shaded regions. Thus, differences among different re- 
alizations can influence only the early evolution of nuclear BHs. At 
z < 1 1 the main BH growth mechanism is gas accretion. 

^ Note that at each redshift/timestep of the simulation, each BH and its host 
halo are the progenitors of the final SMBH and its host galaxy which, by 
construction, will form at redshift z = 6.4. Therefore, at each redshift the 
sum of all single-halo BH masses, as well as of the other quantities, such as 
the DM halo mass, the gas, stars, metals and dust masses, are representative 
of the properties of the final halo at early times (z > 6.4). 
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Figure 4. The black hole accretion rate, Mbh> and black hole mass Mbh 
evolution as a function of redshift (upper and lower panels, respectively), in 
models Q (left panels) and B (right panels). The dot-dashed curves in the 
bottom panels represent the contribution to the BH mass by BH mergers. 




redshift 

Figure 5. The black hole accretion rate, Mbh> and black hole mass, A/bh, 
evolution as a function of redshift in model Q2 for different seed BH masses. 
The magenta (dotted), red (solid) and blue (dashed) lines show the results 
for a single merger tree realization assuming seed BH masses of lO'/i^'Mo, 
10'*/i"'Mo, and lO^ft^'Mo, respectively. 



In Fig. [5] we explore, for a single merger histoiy of model 
Q2, the dependence of the evolution of the BH mass on the as- 
sumed seed BH. We find that models with Mseed ^ 10^ cannot 
reproduce the observed BH mass without increasing the accretion 
efficiency (a). On the other hand, for seed BHs with mass larger 
than ~ lO^'Mg, BH accretion rate and BH mass at the final redshift 
z = 6.4 do not depend much on the adopted seed. This is due to the 
interplay between Eddington limited BH accretion and AGN feed- 
back processes which efficiently regulate BH growth. These results 
justify the seed BH mass adopted in the models (see section [3^.3.1l 



Figure 6. The star formation histories of the host galaxy of the quasar J 1 148 
in quiescent models, Ql, Q2, Q3 (upper panels, from left to right) and in 
bursted models, Bl, B2, B3 (lower panels, from left to right). The an'ow in- 
dicates the upper limit SFR infeiTed from the observation of J1148. The 
SFHs shown in this figure are obtained averaging over 50 random halo 
merger tree realizations and are shown together with the Icr error repre- 
sented by the shaded areas. 



6.2 Star formation histories 

The different SFHs produced by quiescent (Ql, Q2 and Q3) and 
bursted (Bl, B2 and B3) SFH models are shown in Fig. [6] (upper 
and lower panels, respectively), with increasing global star foma- 
tion efficiencies, from low- to high-f,. 

As it can be seen from the figure all SFHs peak at redshift ~ 8, 
with peak amplitudes increasing with /, . The subsequent decrease 
is due to the effect of AGN feedback triggered by a large enough gas 
accretion rate onto the central BH. The galactic scale wind expels 
gas from the halo (see Fig. |7}, reducing the SF rate by a factor of 
~ 2 - 10. The suppression of the SF rate is stronger in Q-models 
where the rate of star formation only depends on the available mass 
of gas at each redshift and does not depend on the merging halos 
mass ratio and redshift distribution. Indeed, such models are only 
mildly affected by the variation among the different merger tree 
realizations, as also demonstrated by the narrow Icr dispersion. 

Conversely, in the B-models, the available gas mass at each 
redshift and thus the shape of the SFH, is determined by the en- 
hanced starburst efficiency which depends on the parameter cTburst- 
The number, intensity and redshift distribution of the starbursts 
strongly reflect the hierarchical merger history of the hosting DM 
halo. In fact, for the same set of parameters, each model shows a 
wide dispersion among the results obtained by individual merger 
histories, with their specific redshift distribution of major mergers 
which trigger the starbursts. The individual SFHs show very differ- 
ent features, such as the onset of starbursts or the redshift at which 
the last major merger occurs. This is reflected in the large 1 cr dis- 
persions which characterize B-models. 

We note that, only model B2 and high-/, models predict av- 
erage SFRs at redshift z = 6.4 in the range (200 - 900) Mg/yi, in 
agreement with what inferred from observations of J1148, as dis- 
cussed in section|2] 
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6.3 Chemical evolution of the host galaxy 

The evolution of the total mass of gas, stars, metals and dust 
(obtained, at each redshift, summing over all existing progenitors) 
for the host galaxy of J 11 48 predicted by the six models is 
shown in Fig. [T] We can immediately see that, independently of 
the assumed SFH, the mass of stars monotonically increases at 
decreasing redshift. On the contrary, the mass of gas and metals 
reach a maximum value at redshift z ~ 8, and then rapidly decrease 
as a consequence of the strong AGN feedback (see also section 
16.2b . It is interesting to observe that, even for the evolution of 
the physical properties, B -models (right panels) show a higher 
dispersion among different hierarchical merger histories than Q- 
models (left panels). We can furthermore note that, except for the 
gas mass which is reproduced by construction, the final properties 
of the host galaxy strongly depend on the assumed SF efficiency, /, . 

Low-f, models are the only ones which predict a final stellar 
mass in good agreement with that estimated from the observed 
dynamical and molecular gas masses of J1148 (M, ~ 4 x 10'° 
Mq). However, both fail to reproduce the final mass of dust and 
metals, which are under-predicted by a factor of 2 - 4. This result 
implies that, independently of the SFH, models which reproduce 
the observed stellar mass are not consistent with the observed 
chemical properties, at least if an ordinary IMF is assumed (see 
section [6!4t . 

Intermediate- f, models predict a final stellar mass (M, ~ 10" Mq) 
larger than observed but still inconsistent with the value expected 
from the local Mbh - M, relation (see Fig. [3}. As a consequence of 
the larger stellar mass, the total mass of metals produced increases 
and it is in (marginal) agreement with the observed value. The final 
dust mass is only reproduced by model Q2. 

High-f, models predict a final stellar mass in agreement with what 
expected from the local Mbh - M, relation [M, ~ (4 - 5) x 10" 
Mq] and reproduce the observed mass of metals. Note that, despite 
the large amount of metals available to accrete onto dust grains, 
the quiescent SFH model, Q3, does not reproduce the observed 
mass of dust. This result, which is opposite to what found for the 
intermediate-/, case, implies that an increase in the total stellar 
mass by a factor of 4 - 5 does not necessarily produce a similar en- 
hancement in the total amount of dust. The origin of this behaviour 
will be discussed in the next section. 



6.3.1 Evolution of dust 

In this section we investigate the evolution of the three components 
that contribute to the total dust content: (i) the mass of dust grown 
in molecular clouds from seed grains produced by stellar sources 
{MC-grown dust); (ii) the mass of dust produced by AGB stars 
(AGB dust); (Hi) the mass of dust formed in SN ejecta (5A' dust). 
The evolution of each component along with the total mass of dust 
(solid red line) are shown in Fig.[8]for all the models. 

It is clear from the figure that dust grown in MCs from seed 
grains produced by stellar sources is the dominant dust component 
at redshift < 10. Indeed, the mass of dust contributed by stellar 
sources (green dot-dashed lines) is always sub-dominant at redshift 
z < 8 - 10, and thus cannot reproduce the observed dust mass. 

The evolution of the different components depends on the 
SFHs (bursted/quiescent) and star formation efficiencies (low- 
/intermediate-/high-/, models). Focusing on stellar dust only, it is 



low-f^ intermediate-f. high-f^ 
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Figure 8. The evolution of individual dust components as a function of 
redshift for the low-f, Ql, Bl (left upper and lower panels, respectively), 
intermediate- f Q2, B2 (middle upper and lower panels, respectively) and 
high-f Q3, B3 (right upper and lower panels, respectively) models. In all 
panels, the total dust mass (red lines) is compared with (i) the dust mass 
computed without the MC-grown dust contribution (green lines), ( ('(') the SN 
dust (cyan lines) and (Hi) the AGB dust components (pink lines). All lines 
are computed averaging over 50 independent merger trees and the shaded 
regions represent the Icr dispersions for SN and AGB-dust components. 

interesting to note that a factor ~ 10 increase in the final mass of 
stars (going from low- to high-/, models) does not imply a simi- 
lar increase in the final dust mass: the average dust mass produced 
by stellar sources (green lines) increases less than a factor 3 and it 
is almost independent of the SFH. Conversely, the relative impor- 
tance of AGB stars and SN depends both on the SFH and on the star 
formation efficiency. In model B2 and high-/, models, AGB stars 
dominate the stellar dust production in the last [150 - 200] Myr of 
evolution, after the peak of star formation and when the effects of 
SN and AGN feedback lead to a strong reduction of the available 
gas mass. Interestingly, the mass of dust produced by AGB stars 
is found to increase in models with a higher star formation effi- 
ciency while a similar effect is not observed for SN dust. Indeed, at 
redshift z = 6.4 the mass of SN dust in high-/, models is compara- 
ble or smaller than in low- and intermediate-/, models. This effect 
can be clarified by the following very simple analytic argument: as- 
suming that SN are the only dust producers, we can integrate eg. 1181 
neglecting the last two terms and applying the instantaneous recy- 
cling approximation to obtain the "effective" dust yield, 

yfsN = ^^^^ ~ Jd,SN - Zd(l + qM„,^,„/snA/'sn) (19) 

where y^sN is the IMF-integrated yield of all stars with masses 
> 8Mq, A/sn is the number of SN formed per unit mass of stars and 
the other symbols have been introduced in section [3.3.41 Thus, as 
the dust-to-gas ratio, Z^, grows in time, the effective dust yield de- 
creases for the combined effect of astration and destruction due to 
SN shocks. The larger star formation efficiencies in high-/, models 
lead to higher dust-to-gas ratios (rapid early dust enrichment and 
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Figure 7. The ISM chemical evolution for the quasar J1148 host galaxy in different models. The lines show the redshift evolution of the average mass of gas, 
stars, metals and dust. Upper panels: low-/» models Ql (left) and Bl (right); middle panels: intermediate-/, models Q2 (left) and B2 (right); lower panels: 
high-/, models Q3 (left) and B3 (right). In all panels solid lines represent the average over 50 merger tree realizations and the shaded regions are the Icr 
dispersion. The arrow is the lower hmit to the gas mass inferred by observations. The circle indicates the final stellar mass with eiTor bars accounting for the 
uncertainties on the estimates of the dynamical mass. The triangle shows the expected average dust mass with eiTorbars accounting for the range of values 
given in the literature (see Table[T). Finally, the filled rectangle represents the estimated metallicity (see text). 
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Figure 9. Time evolution of the accretion (dashed lines) and destraction 
(solid lines) timescales for quiescent models, Ql, Q2, Q3 (upper panels, 
from left to right) and bursted models, Bl, B2, B3 (lower panels, from left 
to right). 




mo'del qr 
. PopIII ■ PopII/L 
PopIII 




Figure 10. The contribution of the first stai's to chemical enrichment of the 
host galaxy ISM as predicted by model Ql. Left panel: the SFH of Pop 
III stars (dotted line) as a function of redshift is compared with the total 
one (solid line). Right panel: the mass of metals (green dotted line) and dust 
(magenta dotted line) produced by Pop III stars at each redshift is compared 
to the total ones (solid ones). The shaded areas represent the corresponding 
Icr dispersion of each quantity. 



in models B2, Q3 and B3 it is limited by destruction. 



higher gas consumption) and are both characterized by lower effec- 
tive dust yields with respect to models with lower /, . For a typical 
dust yield adopted in the present study (see Fig.ll lb. it is found that 
yTsN ~ ^ when Z^ ci ~ ■ Thus, an order-of-magnitude estimate 
of the maximum mass of dust produced by SN can be obtained as 
Zrf.cr X Mgas.max ~ 5 X 10^ Mq where we have assumed the maximum 
mass of gas shown in Fig.|7] This upper limit is consistent with the 
evolution shown in Fig. [8] 

Since AGB stars evolve on longer timescales, their net con- 
tribution is less affected by astration/destruction and the final mass 
of dust they produce increases with the increasing stellar mass (star 
formation efficiency). However, this can only partly compensate for 
the reduced SN contribution, resulting in total stellar dust masses 
which vary less than 30% among different models and never ex- 
ceed 3 X IO'Mq. We conclude that in all models MC-grown dust 
is required in order to reproduce the observed dust mass. It was al- 
ready suggested, although not confirmed, that dust grown in MCs 
could give an important contribution to the observed dust mass in 
Jl 148, given the large molecular gas mass observed in this systems 
(Draine 2009). We find that accretion of dust in MCs can easily 
increase the stellar dust mass by a factor ~ 10 or more, reaching 
final values up to ~ 2 x IO^Mq, the minimum required to match the 
observations. 

Whether or not a model can reach this limit depends on the 
relative importance of the timescales for dust destruction and ac- 
cretion, Tj and Tacc- In Fig.|9]we show these two timescales as a 
function of time for all the models discussed so far: in all cases 
is Tacc « Teh with the exception of the latest stages of the evolu- 
tion when, despite the decline of the supernova rate, gas and metal 
ejection due to AGN feedback conspire to decrease the destruction 
timescale (less supemovae are required to shock the remaining gas) 
and increase the accretion timescale (less metals are available to be 
accreted onto dust grains, hence /^ep ^ 1 in ea.ll2t. 

As expected, is progressively shorter in models where the 
increase of the star formation efficiency leads to larger supernova 
rates. Thus, in models Ql, Bl, and Q2 the final dust mass is 
controlled by accretion (limited by the available mass of metals). 



6.3.2 The contribution of Population III stars to metal and dust 
enrichment 

In this section we analyze the contribution of Pop III stars to metal 
and dust enrichment of the ISM of the QSO host galaxy. We discuss 
the results obtained in model Ql since other models provide very 
similar conclusions. 

The Pop III star formation rate as a function of redshift is 
shown in the left panel of Fig. [TO] (red dotted line with I - cr dis- 
persion). The solid line represents the total (Pop III and Pop II/I) 
star formation rate (same as in the left upper panel of Fig.|6]. As it 
can be seen from the figure, the formation of Pop III stars is dom- 
inant only in the earliest phases of the evolution (z > 16.5) when 
the gas metallicity in progenitor star-forming halos is sub-critical 
(Z < Zct, see section [T2t . In these halos, the gas is promptly en- 
riched by PISNe to metallicities Z > Zc triggering the formation 
of lower mass stars; the formation of Pop III stars is still active 
down to redshift ~ 9.8. However, such starbursts typically occur in 
a small number of halos at each redshift and the resulting Pop III 
star formation rate is orders of magnitude less efficient with respect 
to the total one. 

The right panel of Fig. [TOl shows the mass of metals and dust 
contributed by Pop III stars at each redshift (green and magenta dot- 
ted lines, respectively) with shaded areas indicating the correspond- 
ing 1(T dispersion. PISNe provide up to 100% of the total mass of 
metals and dust at z > 16 but their contribution rapidly drops with 
time and at the final redshift less then [5 - 10]% of metal/dust mass 
is released by Pop III stars. 

Thus, given their negligible contribution, we conclude that the 
choice of the characteristic mass for these stars does not influence 
the chemical evolution of the QSO host galaxy. This result is not 
surprising since, as already emphasized by Schneider, Venkatesan 
& Ferrara (2004), QSOs form in overdense regions of the Universe 
where the Pop Ill/Pop II transition occurs promptly and the enrich- 
ment is dominated by Pop II/I stars. 
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6.4 Chemical evolution with a Top-heavy IMF 



We have shown that low-/, models, which reproduce the final stel- 
lar mass inferred from observations of Jl 148, fail to reproduce both 
the observed dust and metal masses (see Fig.|7] upper panels). De- 
pending on the adopted SFH, a factor of ~ 3 (Q2 model, Fig. [7] 
middle left panel) to ~ 10 (B3 model, Fig.|7] lower right panel) 
higher final stellar mass is required to reproduce the observed 
chemical properties of Jl 148, if a Larson IMF with md, = 0.35 Mq 
is adopted. 

In this section, we investigate how these results depend on 
the stellar IMF. Fig. [TT] shows the evolution of total metal (up- 
per curves) and dust (lower curves) yields as a function of stel- 
lar age, assuming that all stars are formed in a single starburst 
at ? = with solar metallicity (see eq. 1 in V09). Yields are 
computed for a Larson IMF with different characteristic masses, 
rrich = 0.35, 1,5, IOMq, and for a. flat IMF, (f>(m) oc ra"'. This latter 
IMF has been proposed to reconcile predictions of semi-analytical 
models with the observed number counts of submillimeter galax- 
ies at high redshift (e.g. Baugh et al. 2005, but see also Hayward 
et al. 2010). IMFs are all normalized to unity in the mass range 
[0. 1 - 100]Mo and the vertical line indicates the age of Universe at 
redshift z = 6.4. 

We find that a Larson IMF with ra^/, = SMq provides the most 
favourable conditions for larger amount of dust and metals to be 
contemporarily produced by a single stellar population: a factor of 
~ 3 larger metal yield is produced when nid, is increased from our 
standard value 0.35Me to m^, = (5 - 10)Mo. The highest amount 
of dust per unit stellar mass formed (~ 10"^) at z = 6.4 is obtained 
when a Larson IMF with nici, = 5Mo is adopted, which still in- 
cludes a non negligible contribution from AGB stars, comparable 
to that of SNe (see also V09). 

Fig. [12] shows the resulting chemical properties for low-/, 
models when a Larson IMF with md, = 5Mq is adopted. With this 
top-heavy IMF, a factor of ~ 2 - 3 larger dust and metal masses 
at the final redshift z = 6.4 are produced compared to the same 
models with a standard IMF with rtich = 0.35Mo (Fig. [7] upper 
panels). A higher characteristic mass, m^,, implies a larger num- 
ber of intermediate- and high-mass stars. In fact, supernova rates 
comparable to those found in intermediate-/, models, which are 
characterized by a higher star formation efficiency, are produced. 

Since SNe have the dual role of both producing and destroy- 
ing dust, the higher dust yield leads only to a moderate increase in 
the total mass of stellar dust. In Fig.[T3] we show the time evolution 
of individual dust components; the change in the IMF causes SNe 
to dominate stellar dust production at all redshifts. Yet, the total 
mass of stellar dust is still ~ 2 x 10' Mq, consistent with the argu- 
menj^ given by eq.[T9l The larger mass of metals available allows 
~ [2 - 3] X 10^ Mo mass of dust to grow in MCs; by the end of the 
simulation, an equilibrium in reached where ~ t^^^ ; the onset of 
this equilibrium limits the mass of dust. 

Thus, the observed chemical properties of J1148 seem to re- 
quire either models with a factor 3-10 larger stellar mass, which 
would shift this high redshift QSO close or onto the observed lo- 
cal stellar bulge-black hole mass relation, or models with a non- 
standard IMF, in which stars form according to a Larson IMF with 
a high characteristic mass, = 5Me. 



It is important to note that although the SN dust yield ydSN is higher 
when nich = SMq, the number of SN formed per unit mass of stars, A^sN, is 
also higher resulting in a Zj ci- comparable to the value found for a standard 
IMF. 




0.1 
Age (Gyr) 



Figure 11. Evolution of the metal (upper Hnes) and dust (lower lines) yield 
as a function of stellar age. All stars are assumed to form in a single burst 
at f = with solar metalHcity with a Larson IMF with characteristic masses 
nici, = 0.35, 1,5, IOMq and with a flat IMF (see text). The IMFs are normal- 
ized to unity in the mass range [0.1 - lOOlM©. The vertical line indicates 
the age of Universe at redshift z = 6.4. 





Figure 13. The same as in Fig.[8]but for low-/, models with a top-heavy 
IMF (Larson IMF with nici, = 5Mo). 

7 DISCUSSION 

Several theoretical models aimed to investigate the formation and 
evolution of dust in high redshift (z > 5) galaxies and quasars have 
been developed (e.g. Morgan & Edmunds 2003; Hirashita & Fer- 
rara 2002; Dwek et al 2007; V09; Gall et al. 2010; Dwek & Cher- 
chneff' 2010; Gall et al. 201 1; Mattson 201 1). Models which repro- 
duce the observed dust mass in J 1148 using SNe as the primary 
sources require high dust yields (0.3Mo - IMq), special SFHs (all 
stars formed in a starbursts with age 0.3 - 100 Myr), top-heavy IMF 
and/or moderate dust destruction by interstellar shocks (Dwek, Gal- 
Uano & Jones 2007; Dwek & Cherchnefi^ 2011; Gall et al. 2010, 
2011). The high dust yields required by these theoretical calcu- 
lations are not supported by observations of local SNe (Wooden 
et al. 1993; Elmhamdi et al. 2003; Pozzo et al. 2004; Ercolano et 
al. 2007; Meikle et al. 2007; Kotak et al. 2009) and SN remnants. 
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Figure 12. The same as in Fig|2]but for low-f, models with a Larson IMF with niri, = 5Mq. 



Among the latter class, the best studied source is Cas A which has 
been recently observed with Spitzer (Rho et al. 2008) and Herschel 
(Barlow et al. 2010), yielding an estimated dust mass in the range 
2 X IO^^Mq < Mjust < 7.5 X IO^^Mq (see however Dunne et al. 
2009, who proposed the presence of > 0. IM© of dust in Cas A 
based on submm polarization observations with SCUBA). 

Additional dust sources have been also proposed. In V09 we 
pointed out that AGB stars give an important contribution to the ob- 
served dust mass at z > 6 using a SFH for the host galaxy of Jl 148 
as predicted by numerical simulations (Li et al. 2007). This finding 
has been recently confirmed by Dwek & CherchnefF (201 1) and by 
Pipino et al. (2011); the latter also show that, on galactic scales, 
dust produced by the AGN itself - a mechanism known as smok- 
ing QSOs (Elvis et al. 2002) - is negligible when compared to that 
of stellar sources. Quantitative comparisons among various stud- 
ies are difficult due to different prescriptions for computing dust 
formation and evolution (SFH, IMF, dust and metal yields, dust de- 
struction efficiencies). However, all the investigations published so 
far underline the important role played by the adopted SFH and 
IMF 

In this work we have shown that MC-grown dust is the dom- 
inant component even at early cosmic epochs, z > 6, while stellar 
sources, SNe and AGB stars, can not account for the total amount 
of dust observed in Jl 148. Their relative contribution is determined 
by the assumed SFH and IMF. This result differs from what was 
inferred by the simple toy-model presented in V09 for two main 
reasons: the peculiar input SFH and the close-box approximation 
made in V09. Specifically, the redshift of the last major burst of 
star formation and its intensity appear to provide the optimal condi- 
tions for dust production by stellar sources. In turn, our new results 
show the importance of a statistical analysis, where a sufficiently 
large sample of SFHs (merger tree realizations) can be investigated. 
In addition to that, galactic outflows determine the average ISM 
dust-to-gas ratio and have an indirect effect on the dust destruction 
timescales. 

It has been shown in V09 that higher dust masses can be pro- 
duced decoupling the contribution of AGB stars from that of SNe. 
Indeed, we find that the largest stellar dust masses are found in 
some realizations of model B3 where the epoch and amplitude of 
the last major burst maximize dust production by AGB stars. Yet, 



when the average over many independent realizations is consid- 
ered, the effect is partly washed out (see Appendix A). 

Exploring different possible star formation histories for the 
QSO host galaxy, we find that assuming a Larson IMF with 
nich = 0.35Mo, the observed chemical properties (dust mass and 
gas metallicity) of J 11 48 seem to require a final stellar mass which 
is a factor 3-10 higher than what inferred by observations of the 
dynamical and molecular gas masses. It is well known that obser- 
vations of the dynamics and spatial distribution of molecular gas 
in J1148 (Walter et al. 2003), together with those which apply to 
a larger sample of QSOs with 5 < z < 6 (Wang et al. 2010) seem 
to suggest a BH/stellar mass ratio which is 10-20 times higher than 
the present-day value, implying that the central BHs in these sys- 
tems may assemble before their stellar bulge. Observations of dy- 
namical and molecular gas are limited to the central 2.5 kpc of the 
host galaxy; a factor 10 increase in the estimated stellar mass could 
be made consistent with the observations simply assuming that the 
stellar bulge extends out to 25 kpc for realistic radial density dis- 
tributions. More importantly, the required stellar mass - based on 
chemical properties of the host galaxy - would shift the data point 
of J 1148 on the BH-stellar mass relation much closer to the ob- 
served correlation in the present-day galaxies, at least within the 
observed scatter. 

Alternatively, all the observed properties of Jl 148 and its host 
galaxy can be reproduced if Pop II stars are assumed to form ac- 
cording to a Larson IMF with hj^;, = 5Mq, independently of the 
assumed SFH. Since SNe both form and destroy dust grains, the 
larger dust yield due to higher supernova rate leads only to a mild 
increase in the total dust mass produced by stellar sources. Yet, 
the larger metal yields lead to larger masses of gas-phase metals, 
which allow seed dust grains to efficiently grow in MCs. Indeed, 
physical conditions present in the ISM of QSO host galaxies at 
high-z may favour a shift in the fragmentation mass scale of star 
forming clouds, hence in the characteristic mass of the stellar IMF. 
Schneider & Omukai (2010) have studied the response of the ther- 
mal Jeans mass at cloud fragmentation to the joint effect of metal- 
line cooling, dust-cooling and the cosmic microwave background 
(CMB) in the ISM of high-redshift galaxies. They show that when 
z > 5 and the metallicity is Z > IO^^Zq, heating of dust grains by 
the CMB shifts the Jeans scale to larger masses. If no dust grains are 
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present in star forming cloud, the effect is even more dramatic and 
only massive and very massive stars are predicted to form (Smith 
et al. 2009; Jappsen et al. 2009). In addition to the effect of the 
CMB, MCs in circumnuclear starbnrst regions may be character- 
ized by gas temperatures ~ 100 K and densities of a few 10^ cm"'; 
under these conditions numerical simulations predict that the re- 
sulting mass spectrum of gravitational condensations (a proxy to 
the IMF) is top-heavy, compared to that in the solar neighbourhood 
(Klessen, Spaans & Jappsen 2007). 

In the present study (as in VOQ) SN dust yields are obtained 
using one of the models proposed by Bianchi & Schneider (2007) 
where only 7% of the dust formed in SN ejecta survive the passage 
of the reverse shock. According to this model, ~ [0.5-5]x 10"^) Mq 
of dust per SN is injected in the ISM (depending on the progenitor 
star mass and metallicity), in good agreement with dust observed 
in young SN remnants. It is important to discuss how a different 
choice would affect the results present in previous sections. Indeed, 
if we adopt an alternative model proposed by Bianchi & Schneider 
(2007), where the fraction of SN-condensed dust which is injected 
in the ISM is 20%, we find that the final dust masses are similar- 
ily shifted upward by a factor ~ 3, bringing models Ql and BI in 
agreement with the observational data point. However both models 
severely under-predict the mass of metals. Thus, although a higher 
SN dust yield leads to a larger total dust mass (closer to the ob- 
served value), a larger final stellar mass or a top-heavy IMF are 
still required to reproduce all the observed properties of J 1 148. 

A final remark concerns the adopted prescription for the 
AGN feedback. Feruglio et al. (2010) have recently detected the 
broad wings of the CO line in Mrk 231, the closest known QSO, 
which has enabled to trace a giant molecular outflow of about 
700 Mg/yr; the inferred kinetic energy in the molecular outflow 
is 1.2 10'*'* erg/s, corresponding to a few percent of the AGN bolo- 
metric luminosity (5 10*^ erg/s). This is consistent with what we 
find in our simulated models. In fact, the formulation adopted for 
the BH accretion rate ensures that BH feedback is active over al- 
most the entire lifetime of the host halo (since the computed rate is 
almost always above the critical threshold), resulting in an efficient 
gas ejection at lower redshift (z < 7 - 8) when the BH mass and 
its gas accretion rate are higher. This, in turn, is regulated by the 
choice of the parameters a and ew,AGN- However, it is important to 
note that the relative role of starburst and AGN activity in regulat- 
ing/quenching star formation and BH-growth, as well as the mech- 
anisms (energy/momentum deposition) driving these feedback pro- 
cesses, are still largely debated (see e.g. Bhur et al. 2010; Hopkins 
201 1; Hopkins, Quataert & Murray 201 1 for a different view). 



8 SUMMARY AND CONCLUSIONS 

Here we have presented a semi-analytical model which allows to 
(ij model the evolution of a QSO host galaxy and its central su- 
permassive black hole, predicting a range of possible SFHs and 
to (/(') explore the relative importance of stellar sources of dust 
(SN and AGB stars) and of dust grown in MCs at these early cos- 
mic epochs. GAMETE/QSOdust follows the build-up of the cen- 
tral SMBH through mergers and gas accretion, starting from seed 
black holes at very high redshifts. The corresponding assembly of 
the host galaxy is followed in a self-consistent way and the co- 
evolution of these two fundamental components is controlled by 
AGN feedback in the form of a galactic-scale wind, assuming that 
0.5% of the rest mass energy of the gas accreted onto the BH is 
thermally coupled to the ISM. We have performed 50 random hi- 



erarchical merger histories of a lO'^M© dark matter halo, which is 
believed to host a > IO'Mq SMBH, through a binary Monte Carlo 
algorithm with mass accretion based on the EPS theory. 

Metal and dust enrichment is followed evolving all progeni- 
tor galaxies on the corresponding stellar lifetimes, considering the 
contribution from Pop III (PISNe) and Pop II stars (AGB and SNe), 
and the subsequent dust grains processing (destruction and growth) 
in the ISM. Metal and dust yields are taken from theoretical mod- 
els for stellar nucleosynthesis and dust nucleation in SN ejecta and 
AGB stellar atmospheres. Different SFHs for the QSO host galaxy 
have been explored: in quiescent models, the efficiency of star for- 
mation is independent of the galaxy mergers mass ratio, while in 
bursted models the star formation efficiency is enhanced during ma- 
jor mergers among progenitor galaxies. The free parameters of the 
models, namely the efficiency of star formation, BH accretion and 
AGN-driven wind have been chosen to reproduce the BH and gas 
mass inferred from observations of Jl 148 at redshift z = 6.4. 

The main results of this study can be summarized as follows: 

• the chemical properties of the host galaxies of high-z QSOs 
are a powerful probe of their past evolution. The dust mass pro- 
duced by stars depends on the shape of the star formation history, 
which is controlled by the evolution of the gas content through 
mergers and outflows driven by SNe and by the AGN. The total 
mass of metals is a tracer of the integrated star formation history 
(total stellar mass) and it is sensitive to the strength and frequen- 
cies of metal-rich outflows. 

• Population III stars give a negligible contribution to chemical 
enrichment of the ISM. In fact, the large metal and dust yields ex- 
pected from theoretical models of PISN lead to rapid enrichment 
of the ISM to metallicities Z > = [10~^ - 10"'']Zo in progenitor 
galaxies of the final galaxy host, enabling the formation of Popula- 
tion II/I stars. 

• If Population II/I stars form with a standard IMF and a char- 
acteristic stellar mass of rrich = 0.35 Mq, as in the present-day Uni- 
verse, the final masses of metals and dust are lower than observed. 
Depending on the adopted SFH, a final stellar mass which is a fac- 
tor of 3 - 10 larger than the value inferred from observations of 
the dynamics and spatial distribution of the molecular gas content 
of Jl 148 (Walter et al. 2003) is required. The required stellar mass 
would shift the position of J 1148 onto the locally observed corre- 
lation between the SMBH and stellar bulge masses, at least within 
the observed scatter. 

• If Population II/I stars form with a top-heavy IMF and a char- 
acteristic stellar mass of nich = 5 Mq then the observed chemi- 
cal properties can be reconciled with the inferred stellar mass. The 
physical conditions of star forming regions in these high-z systems 
may in fact suppress fragmentation, due to higher temperatures in 
MCs and heating of dust grains, favouring the formation of larger 
stellar masses. 

• A statistical analysis of the SFHs, with many independent hi- 
erarchical merger trees of the dark matter halo hosting J1148, has 
allowed us to conclude that although SNe dominate the early dust 
enrichment, AGB stars contribute at z < 8 - 10; the latter contribu- 
tion depends on the shape of the SFH and on the adopted IMF for 
Population II/I stars. 

• We find that stellar sources produce a total dust mass which 
never exceeds ~ 2 - 3 x IO^Mq; in fact, unless the redshift and in- 
tensity of the last major burst of star formation allows to efficiently 
decouple the contribution of SNe and AGB stars, the effective dust 
yield is progressively reduced due to the joint effect of astration 
and dust destruction in the ISM by SN shocks. This conclusion is 
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model 






^burst,max 


a 




B4 


0.1 


0.1 


4.0 


233 


5 X 10"^ 


B5 


0.1 


0.25 


1.6 


250 


5 X 10-3 



Table Al. Model parameters of the two additional bursted SFH models, B4 
and B5. Column descriptions are the same as in Table |2] 

independent of the adopted stellar IMF but depend on the adopted 
stellar dust yields. 

• In all the models we have investigated, the final dust mass can 
be reconciled with the observed value, Mj^s, = [2 - 6] x 10^ Mg 
via grain growth in molecular clouds which appear to be important 
even at these high redshifts, consistent with the large molecular gas 
mass observed in JI 148. The final mass of dust is thus controlled by 
two fundamental timescales: the timescale for grain growth (which 
become progressively longer due to rapid depletion of metals in 
dust grains) and the timescale for dust destruction by SN shocks 
(which depends on the SFH through both gas consumption/ejection 
and the SN rate). 
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Figure Al. SFHs in high-/» models B3, B4 and B5 for a single merger tree 
realization (thin lines). The panels show the eft'ects of dift'erent values of 
the parameter CTi,urst'- 0.05 (B3, left panel), 0.1 (B4, central panel) and 0.25 
(B5, right panel). Thick solid lines refer to the average over 50 merger tree 
realizations (see text). 
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APPENDIX A: DEPENDENCE ON THE STARBURST 
FREQUENCY 

With the aim of investigating the dependence of our results on the 
number and intensity of starbursts, we have explored two additional 
high-/, SFHs, models B4 and B5 (see Table lAlb . These models pre- 
dict a total stellar mass comparable to model B3 but with increasing 
values of CTburst which imply SFHs characterized by an increasing 
number of (less efficient) bursts of star formation. The resulting av- 
erage SFHs are shown in Fig. lAll (thick lines and corresponding 
l-o" dispersion). The two SFH models have progressively less in- 
tense but more frequent bursts with respect to model B3. These dif- 
ferences are particularly evident when we compare the predictions 
of a single merger tree realization; the thin lines in Fig. lAll show 
the predicted SFHs in models B3, B4, and B5 for the same merger 
tree realization. Comparing with the corresponding thick lines in 
each panel, which represent the average SFH over 50 merger tree 
realizations, it becomes clear how the averaging procedure partly 
cancel the redshift-dependent effects. 

The quasar host chemical evolution predicted by the additional 
models presented here, B4 and B5, are shown in Fig. IA2I Even if 
the evolution at the highest redshifts is (slightly) different (the pre- 
dicted masses are slightly higher and the curves are progressively 
smoother), the final masses of gas, stars and metals are very similar 
to what predicted by model B3 (right lower panel in Fig.|7). On the 
other hand, we find that model B4 is in agreement with the observed 
dust mass within the Icr dispersion on the average evolution (i.e. at 
least in one of the 50 merger histories the associated SFH is able 

' http://wiki.arcetri.astro.it/bin/view/DAVID/WebHome 



Figure A2. Evolution of the chemical properties (gas, stars, metals and 
dust) for bursted SF models B4 and B5. Line colours and observational 
constraints are the same as in Fig. [7] 

to reproduce the observed properties of the quasar host), whereas 
model B5, which is characterized by a SFH with less intense and 
more frequent bursts, is not able to reproduce the observed dust 
mass. 

This analysis shows that: 

(i) By averaging the results over a set of merger tree realizations 
the specific features of each single SFH, such as the number, in- 
tensity and redshift distribution of starbursts and the intensity and 
position of the burst triggered by the last major merger, are partly 
washed out. 

(ii) An increasing number of starbursts does not produce a larger 
final dust mass. The higher is the number of bursts (over the whole 
host galaxy evolution), the smaller is their intensity, approaching a 
"quiescent like" SFH (see model B5); the final dust mass is strongly 
limited by destruction (see section l63Tt . 
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